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Abstract 

We describe a system composed of two coupled optical cavity modes with a coupling modulated 
by a bulk mechanical resonator. In addition, one of the cavity modes is irreversibly coupled to a 
single photon source. Our scheme is an opto-mechanical realisation of the Jaynes-Cummings model 
where the qubit is a dual rail optical qubit while the bosonic degree of freedom is a matter degree of 
freedom realised as the bulk mechanical excitation. We show the possibility of engineering phonon 
number states of the mechanical oscillator in such a system by computing the conditional state of 
the mechanics after successive photon counting measurements. 

PACS numbers: 42.50.Wk, 42.50.Lc,42.50.Dv 



* Author to whom any correspondence should be addressed. Electronic mail: s.basiri@uq.edu.au 



1 



I. INTRODUCTION 



Quantum opto-mechanics [U-S] provides an exciting context in which to investigate the 
interaction of hght and the collective motion of bulk mechanical systems in the quantum 
regime. More generally, quantum opto-mechanics provides a good example of an engineered 
quantum system: a meso/macroscopic device specifically engineered so that collective de- 
grees of freedom can be subject to coherent quantum control and measurement. Recent 
achievement of the quantum ground state of a mechanical resonator in electro [5l E] and 
opto-mechanical systems [7] has meant that we can now begin to observe and eventually 
manipulate quantum states [H [9] at the mesoscopic scale. A recent opto-mechanical exper- 
iment towards this direction has been carried out illustrating an efficient quantum interface 
between optical photons and mechanical phonons [TDj . An obvious next step is understand- 
ing and implementing single photon optical non linearities in such systems. Recent results 
for single photon opto-mechanics describe novel features such as the photon blockade ef- 
fect p!T| [T2] and cavity resonance shift [13] , preparation of non gaussian mechanical states 
[m [15] , and mechanical superpositions [I6] . 

The model of this paper is based on two optical cavity modes with a coupling modulated 
by a mechanical degree of freedom. This can be realised in a number of different ways, but 
we will focus on the scheme of Chang et al. [17]. In this scheme a coherent exchange of 
photons between two photonic crystal defect cavities is modulated by a common mechanical 
degree of freedom. This implementation leads to a large coupling strength between the 
optical and mechanical modes and is a promising direction to achieve strong coupling at the 
single photon level. Another implementation could be based on a single bulk flexural mode 
driven by the opposing radiation pressure forces of two optical cavity modes. If the cavity 
modes are coupled, transformation to normal modes leads to a model in which the normal 
mode coupling is modulated by the mechanical displacement jTg| . 

In the present work we consider single photon driving of one unit in the opto-mechanical 
crystal of Chang et al., and show how this system may be mathematically described using 
the Jaynes-Cummings model [19] when the optical modes are excited by single photon states. 
A similar analogy to the Jaynes Cummings Hamiltonian has been studied previously in an 
atom assisted cavity optomechanical system in [20] . In our study the bosonic component of 
the Jaynes-Cummings model becomes the mechanical degree of freedom while the two-level 
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component of the Jaynes-Cummings model is a dual-rail single photon qubit. The Jaynes- 
Cummings analogy can be exploited to realise a measurement of the phonon number of the 
mechanical degree of freedom following a modified version of the measurement scheme of 
Guerlin et al. [21]. Additionally, we condition the system on photon counts and calculate 
conditional states of the mechanical resonator. 

A schematic diagram of the model we treat is given in figure [1} Two optical cavity modes 
are coupled in such a way that the coupling is proportional to a mechanical displacement. 
This leads to a Raman process in which photons are exchanged between the cavities by 
absorbing or emitting a phonon. The electrical field in each cavity is described by a single 
mode with photon annihilation and creation operators ai,a\ {i = 1,2). The mechanical 
degree of freedom is described as single harmonic oscillator with phonon annihilation and 
creation operators 6,6^. We will assume that one optical cavity, 02, is excited by emission 




FIG. 1. A schemata for an opto- mechanical system in which two single photonic cavity modes ai, 02, 
coherently exchange photons at a rate proportional to the displacement of a collective mechanical 
degree of freedom, b. The single photon source cavity is labelled c and excites cavity 02 irreversibly 
through the presence of a circulator, labelled cir. Single photon counters are labelled Di, D2- 

from a single photon source modelled as another optical cavity (the source) prepared with 
exactly one photon. The single photon excitation of 02 is an irreversible process due to 
insertion of a circulator between the source and cavity 02- The injected photon can be 
reflected from cavity 02 or absorbed and then remitted, to be detected at a single photo 
counter D2. We also allow for the possibility of emission from cavity ai which can be 
monitored by detector Di. 
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The Hamiltonian for the opto-mechanical system, comprising the two cavity modes and 
the mechanical degree of freedom, may be written as [T71 [TS] 

H = huJia[ai + hu2a!2a2 + hujmh'^h + hg{h + })^){a\a2 + OiOg), (1) 

where uJi/2 is the resonance frequency of each optical cavity, Um is the mechanical resonance 
frequency and g is the rate of coherent coupling between the mechanics and the two optical 
cavities. Other realisations of such a three mode optomechanical interaction include [22] 
using the membrane in the middle model and also |23] in the context of an opto-acoustic 
parametric amplifier. We now move to an interaction picture and assume that the system 
is so designed that uj2 = oJi + uJm- The opto-mechanical interaction picture Hamiltonian, 
including only resonant terms, is 

Horn = hg{b^a[a2 + baial). (2) 

This represents a kind of coherent Raman process whereby one photon from cavity 02 is 
transferred into cavity Oi simultaneously exciting one phonon in the mechanical degree of 
freedom. 

The analysis is greatly simplified if we assume that, at most, there is one photon in the 
system at any time. In that case, the opto-mechanical interaction can be regarded as an 
interaction between a qubit and a simple harmonic oscillator with the qubit states defined 
as 

|0) = |l)i|0)2, (3) 
|l) = |0)i|l)2, (4) 

where \n)i are photon number eigenstates for cavity This is the dual rail encoding used 
in linear optical quantum computing schemes [211 EH]. On this restricted subspace we can 
define 

a\a2 = \0){l\=a., (5) 
aial = \l){0\=a+, (6) 

where a± are the usual raising and lowering operators for a pseudo-spin system. The opto- 
mechanical interaction Hamiltonian can then be written in terms of the Jaynes-Cummings 
Hamiltonian 

Hom = hg{ba+ + b^a_). (7) 



Each cavity mode is treated as a single sided cavity with photon decay rates given by 
Ki/2 for each cavity mode respectively. The restriction to a single sided cavity is purely for 
simplicity in presenting our argument. It is straightforward to include loss at a second cavity 
mirror. The additional output port, if not monitored, provides an additional decay channel 
for the single photon. Any photon lost through the unmonitored port is never detected and 
thus leads to an increase in the number of trails required for a successful detection at the 
monitored port. 

The single photon source, modelled as a source cavity prepared with a single photon [S], is 
irreversibly coupled into the input of the cavity 02. The emission rate from the source cavity 
is 7. The input single photon states thus constitute a pulse with an exponential temporal 
profile with lifetime 7"^. Other single photon excitations mechanism could be implemented 
including a quantum dot source embedded in the waveguide as in Schwagmann et al. [26] . 

We model this kind of single photon excitation process using the theory of cascaded 
quantum systems [2H] to give a master equation for the dynamics of the OM system plus 
the source cavity 

^ = -^[i/o™.p]-^^[c4-c^«2,p]+I?[J]p+/ti^^[ai]p+7rn(A^+l)^^[&]p+7™iVI?[6t]p, (8) 

were c, are the annihilation and creation operators for photon number states in the source 
cavity, the superoperator V is defined by 

V[A]p = ApA^ - ^{A^Ap + pA^A), (9) 

and the detection operator is J = ^/7C + ^^02 which consists of the sum of two terms: 
the first term shows that the photon can be detected directly from the source before it is 
transmitted to cavity 02 and the second term shows a photo-detection from within cavity 02. 
We have also assumed that the source cavity is on resonance with the receiver cavity, 02. If 
cavity-2 had two output mirrors, we would include an additional decay term in the master 
equation of the form K2^[a2]p where K2 is the loss rate through the second, unmonitored 
output mirror of this cavity. 

The initial state is taken to be such that at t = there is one photon in the source and 
no photons in either cavity ai or cavity 02, while the mechanical system is in an arbitrary 
coherent state \(3)b = e~'^^'^'^ ^'^^^(3'^ / ^/n\\n)]j where l)^h\n)n = n\n)n. The total initial state 
is thus 

|vl>(0)) = |l)e|0)i|0)2|/3),. (10) 
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Given this initial state, we now ask for the conditional state of the mechanics given that 
the photon is detected at D2 at a time ti. In the ideal case we would like the photon to be 
detected with certainty at D2, however in reality it could fail to be counted due to non unit 
quantum efficiency in the detector, lost either through emission out of cavity ai, scattered 
out of the input and output channels, or perhaps not emitted by the source at all. However 
assuming the ideal case, every photon emitted by the source is counted at D2 or Di. If 
we further assume the decay rate of cavity ai is negligible (although we include it in our 
analysis), every photon emitted by the source is counted so that no information is lost and 
each photon detection event at D2 gives a single bit of information. To gain additional 
information about the mechanics we can simply repeat the process, each time preparing 
a single photon in the source and using the conditional mechanical state obtained by the 
previous detection event. 

II. MEASUREMENT OF PHONON NUMBER 

A. Measurement and the Jaynes-Cummings model 

At first sight the interaction Hamiltonian in Eq.Q does not look as if it could realise 
a measurement of the phonon number operator, for the mechanics. However given the 
Jaynes-Cummings representation in Eq.([7]) we can use the results of the Haroche group [29] 
to see how it may be configured so as to yield information on the state of the mechanical 
system. We will first review the simpler case in which interactions are deterministic. 

Let us suppose that we prepare the two-level system in the state |1) and that the inter- 
action with the bosonic system proceeds for a time r at which point we make an arbitrarily 
accurate measurement of the qubit state. The result of this measurement is a single binary 
number x. The resulting conditional (unnormalised) state for the bosonic degree of freedom 
is 




(11) 



where 




(12) 
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where 6 = gr. It is a simple matter to show that 

E{1) = cos{eVW), (13) 
E{0) = -ib^bb'^)-^^'^ sm{eVm). (14) 

The probabihty to obtain the result x is simply the required normalisation of the unnor- 
malised conditional state 

p{x) = = b{m\EHx)Eix)\m)b. (15) 

As i?"''(x)-E'(x) = 1 this probability distribution is normalised. 



As the measurement operators in Eqs. (13) and (14) commute with the mechanical 



phonon number, this model describes a coarse-grained phonon number measurement; coarse- 
grained because there is only one bit of information per trial. It is not however a quantum 
nondemolition measurement (QND) as the interaction Hamiltonian between the probe and 
the mechanical system in Eq. ([s]) does not commute with the number operator. In fact, as 
we started with an interaction that is linear in the mechanical displacement we could not 
expect a QND measurement of the phonon number which requires an interaction that is 
at least quadratic in the displacement [30j. It is possible to configure the system discussed 
in this paper in the strong dispersive regime so that there is an effective coupling to the 
displacement squared thus realising a QND phonon number measurement [18]. 

Suppose the initial state of the bosonic degree of freedom is a coherent state \ip{0)) = \/3)h, 
where 



oo 



n v'n'. 

71=0 



so that the initial number distribution is the Poisson distribution 



Po(n) = e-l'^l^^. (17) 



If the result of the measurement was x = 1, the number distribution for the conditional 
state now becomes 

Pi{n) = [p{l)r^ cos\9V^^^)Po{n). (18) 

In figure [2] we plot this distribution, together with the initial distribution, for various 
values of the parameter 6. The important feature to note is that for particular values of 6, 
in this case 6 = 7r/6, the distribution is very different from the initial poisson distribution. 



Figure [s] shows the probabihty for this outcome s a function of 6, p{l) = {cos{9y/n + 1)) 
with the average taken over the number distribution prior to measurement. Clearly this 
average is bounded by one and, from the figure, we see that it oscillates around 1/2. This 
indicates that the information provided in a single measurement is at most one bit. 
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FIG. 2. The conditional number distribution after one readout with x = 1 for various values of 9. 
In all cases /? = 3. 
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FIG. 3. The probability to obtain the result x = 1 for /3 = 3 versus 9 = gr. 

We now consider what happens if we repeat the measurement using the conditional state 
from the first measurement as the initial state for the bosonic degree of freedom for the 
next measurement, but change the value of 6. In figure |4] we plot the conditional number 
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distribution for fifty five measurements that all gave the result x = 1, but with different 
values of 6. We see that after the sequence of measurements, for appropriate values of 6, 




FIG. 4. Photon number probabilities versus photon number and number of readout, all with x = 1, 
/3 = 3 and for various values of 9. 

the conditional state can approach a number state. 

To get an idea of how many measurements are required to reach a number state, we can 
consider the case for which the value of 6 in each trial is the same. In such trials, the 
final distribution is given by 

PNin) = 1 . . . , 1)]-^ cos''''ieV^^)Poin). (19) 

where the normalisation is simply the probability for this history of measurement results, 
1 . . . , 1) = cos^^(6'-\A?r+T)Po('^)- The cosine factor is a periodic comb-like dis- 

tribution with respect to n. If we choose 6 to align the comb so that a peak is on the 
mean phonon number, we can estimate the width of the distribution P/v(^)- We thus choose 
6y/^ = vr, and expand the cosine factor around this value to find that the distribution can 
be approximated by a Gaussian with variance 



This indicates that the number of measurements we would need to make to approach a 
Fock state, for which W = 1, scales as N oc n^. This suggests that one should prepare the 
mechanical resonator in a state with a small value of n in order to have a reasonable chance 
of getting to a Fock state before mechanical heating takes over. 

In the Haroche experiments [Mj the values oi 9 = gr are kept approximately constant 
from trial to trial by adjusting the velocity distribution of the atoms passing through the 
cavity. In the model of this paper however the interaction times r are stochastic due to the 
probabilistic character of the photon dwell time in the cavity. However it is clear that the 
interaction times are most likely to be of the order of k^"^ so that the condition 6^/fl = vr 
translates to ~ n j^fi. For reasonable values of n, this does not make it too difficult to 
achieve reasonable values of single-photon optomechanical coupling rate. 

B. Conditional mechanical state based on photo-detection 

1. Without mechanical damping 

In the case of the opto-mechanical model, the conditional state of the mechanics depends 
on the time of detection, ti, which is itself a random variable. This is quite different to the 
simple Jaynes-Cummings measurement model considered in the previous section where the 
interaction time r was under the control of the experimenter. However, in the ideal case, if 
a photon is emitted from cavity 02 it will be detected at As the emission rate is simply 
proportional to the probability of there being one photon in this cavity, a detection at some 
time ti is equivalent to a direct readout of the qubit state |1) in the Jaynes-Cummings 
model with an interaction time ti although, unlike that model, the qubit is destroyed in the 
process. 

However there is an additional new feature. A photon detection at D2 can occur in two 
indistinguishable ways: the photon can be reflected from cavity 02 directly into the detector 
or it can be transmitted from within the cavity after first being absorbed. This feature is 
refiected in the jump operator J as the sum of two terms, J = ^/^c + y/K2(^2- This can lead 
to an interference term in the detection rate. 

We also need to include the possibility that a photon is lost from cavity ai. As there 
is at most one photon in the entire system at any time, a photon lost from ai means that 
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no photon can be counted at D2- Thus loss from the other cavity appears as a non unit 
efficiency in the detection process. One could, of course, insert a detector at the output 
from Oi to herald such erroneous events. If we do not monitor this channel, we will need 
to define a cut-off time: non-detection at D2 up to the cut-off time indicates failure and we 
simply discard that run and start again. 

We ffist compute the conditional state conditioned on no detections, either at Di or 
D21 up to time t. The quantum theory of continuous photon counting [31] shows that this 
(unnormalised) conditional state is determined by 

p(°'°)(t) = 5(t)p(0), (21) 

where the superscript is defined by (ni, 722) where is the count number recorded at detector 
Di. where iS(t)p(0) = p^^'^\t) is the conditional state given no counts up to time t and is 
given by solving 

f^=-r{K~p-~pK^), (22) 
where the non-Hermitian operator K is given by 

K = g{b^a\a2 + baial) — i^f^~{K2{c(i}2 — c^a2)/2 — iJ^ J /2 — iKia\ai/2. (23) 

The normalisation of this state is simply the probability for no counts up to time t, 

p(ni = 0,n2 = 0,t) = tr[p(0'°)(t)]. (24) 

Note, that if the initial state is pure, we need only to solve the effective Schrodinger equation 

- -^Km, (25) 

to give 

p{m(t) = |v^(°(t))(^(°)(t)|. (26) 

We now ask for the conditional state of the system given that one photon is counted at 
D2 in time t to t -|- dt. Such an event means that no photon can have been decayed through 
the output of cavity Oi. This conditional state is 

p(°-^)(t) = J5(t)p(0)jt. (27) 

If the initial state is a pure state, this conditional state is also a pure state 

1^(0.1)) = J|^W(t)) = y^c|V^(°)(t)) + y^a2|V^(°)(t)), (28) 
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FIG. 5. The photon detection rate (sohd hne) as the sum of three terms: source emission rate 
(dotted hne), cavity emission rate (dashed hne), interference emission rate (dashed-dotted hne) 
showing how interference leads to a minimum at finite time. The zero is due to interference 
between the two indistinguishable ways a single photon can be detected: reflected from the input 
mirror or absorbed by the cavity and subsequently re-emitted. The parameters are g = and 
Kl = 0. 

which is superposition of the two ways in which a photon can be counted: direct reflection 
form the cavity or emission from inside the cavity. This leads to an interference term in the 
detection rate 

Roiit) = 7(^(°)(t)|ctc|^(°)(t)) + K2{^p'^'\t)\alaS^'{t)) (29) 
+Vl^2{{i'^'\t)\c^aS^%t))+c.cy 

For example, if Ki = and g = 0, we find that 

where the mean photon number in cavity 02 is 

r^2(t) = ^^^(e-^*/^-e--*/^)^ (31) 

{K2 - 7) 



The first two terms in Eq. ( 30 ) correspond to direct detection from the source cavity and 02 
respectively. The last term represents the interference between these two indistinguishable 
events and leads to a zero in the detection rate at D2 as a function of time (figure 5). 

If the photon decays through cavity ai, it can never be detected at D2. However if we do 
not monitor this output we have no way of knowing when this happens. We thus need to 
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sum over all times at which a photon could be emitted from cavity ai. Note also that once a 
photon is lost the operation S acts trivially as the identity. The (unnormalised) conditional 
state given that one photon was lost from cavity Oi at any time over the interval [0,t) is 

p(i'0)([0,t)) = /ti f dtiaie-'''''p{0)e'''''a\. (32) 
Jo 

The normalisation of this state is the error probability 

Perr{t) = ^1 f dtiti L e-*^*V(0)e^^*^a| 1 , (33) 
Jo ^ ^ 

as this represents the probability that a photon was lost from the system before it could be 
detected at D2. Clearly to keep the error low we need to ensure ni « K2. 

If we expand the initial mechanical state in the eigenstates of 6^6 which we write as \n)b 
the dynamics is closed in the three dimensional subspace spanned by the basis 

{|l)„,|2)„,|3)„} = {|l,0,0,n),|0,0,l,n),|0,l,0,n + l)}, (34) 

where z,n) = \x)c\y)i\z)2\n)b- We can now expand 
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\^^'\t)) = Y,Cm\k)n. (35) 



k=l 



Substituting this into Eq.(25) gives a closed set of equations for the coefficients that can be 



solved analytically. For future purposes we rewrite the initial state in Eq.(lO) in another 
form 

00 

l^(0)) = $^/3n(to)|l)n, (36) 

n=0 

where /3n(to) = Therefore, the initial number distribution is the Poisson distribu- 

tion 

P'nito) = e-'"''/'^- (37) 



As discussed in section II A we need to be in strong opto-mechanical coupling regime 
for which g is of the order of ^2 in order to reach (or get close to) a Fock state within a 
reasonable number of trials. Hence, in units such that ^2 = 1, we have g = 1. We also 
assume that ki = 0.2, 7 = 0.9 and the mechanical damping rate, 7^ << 1^21 so that we can 
neglect the mechanical damping. We consider mechanical damping and thermal effects in 
subsection IIIB 2[ 
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FIG. 6. The photon detection rate for /3 = 2, 51 = 1, K2 = 1 and ki = 0.2. 

Figure 6 shows the detection rate for the first measurement for the given parameters. We 
use the rate function to generate a random detection time ti in a time interval which starts 
just after the minimum and ends close where the detection rate is nearly zero. This choice 
assures us that, with high probability, we are detecting a photon from cavity 02, after it has 
interacted with the mechanical system not one which is reflected of the mirror directly from 



the source. We then substitute ti in Eq.(35) to get the normalised conditional state of the 
system given no counts up to ti 

|,<o,,y> ^ , l^'""-') . (38) 



Applying jump operator J on \i^^^\ti)) we get the conditional state of the system \ip^^\ti)) 

given that one photon is counted in time ti to ti + dti. Normalising this state we get 

00 

\i^^'\h)) = Y,Pn{h)\0)n, (39) 
n=0 

where 

^^^^^^ ^ Pnito)i^cUtl) + y^2cUtl)) ^^Q^ 



Tr=o/3^(^o)lv/7c?(ti) + v^c^(ti 
The phonon number distribution for the conditional state now becomes 

P^Xh) = P^MPin,h), (41) 

where 
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We will repeat the measurement process preparing another single photon in the source 
and using the new initial state Yl'^=o f^n(ti)\^) ■ Given the state of the system after each 
measurement, we can re-calculate the detection rate from which we again sample a random 
detection time. 

After r detection events, the conditional state of the system is 



n=0 



(43) 




FIG. 7. Phonon number distribution function histograms after successive measurements for (3 = 2 
(n = 4), = 1, 7 = 0.9, K2 = Ij /^i = 0.2 and random detection times. As we increase the 
number of measurements, the phonon number distribution evolves from a poissonian distribution 
into number state distributions (a) n=2 and (b) n=3. 

where 



(45) 



The phonon number distribution function after r measurements now becomes 



where 



P(n, tr 



|y7C?(tr) + v^C^(tr)P 



(46) 
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Each measurement provides partial information about phonon number (typically less than 
one but per trial). The procedure can be explained looking at Eq.([45|) and is quite similar 
to the simple model with deterministic interaction times discussed in section II A After 
the r'th detection event, the phonon number distribution is multiplied by a filter function 
P{n,tr) which, for appropriate values of t,., suppresses certain values of n. Continuing 
the measurement process we can get more information leading to a gradual collapse of the 
distribution onto a single number state. Figures 7(a) and 7(b) show the phonon number 
distributions for sixty successive measurements and eighty measurements respectively. The 
number of measurements required to reach a Fock state is at least roughly consistent with the 



prediction of the simple model in II A , i.e. scaling as n^, despite the stochastic fluctuations 
in the interaction times in the optomechanical model. 
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FIG. 8. The evohition of the Wigner distribution function Wr{x,p) with successive measurements 
where r is the measurement number. (3 = 2, g = 1, j = 0.9, K2 = 1, /^i = 0.2 and the same 
detection times as in figure 7(b). Measuring the qubit repeatedly we drive the mechanics from 
a state with a distribution function from a Poisson distribution with mean n = 4 to a number 
eigenstate with exact excitation number of n = 3. 



The figure shows the evolution of the number distribution from Poisson distribution into 
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Fock state at n — 2 and n = 3 for different cfioices of the simulated detection times. Once 
the system collapses to a specific phonon number state, it will remain in that state upon 
further measurements. (Given that we have neglected mechanical dissipation and thermal 
fluctuations). Some simulations do not settle down to a single number state but show 
jumping between two number states so that in some cases we have a competition between 
two n values in number distribution. 

We can illustrate that the steady state conditional states do indeed tend to number states 
by computing the Wigner function. Figure 8 shows the Wigner functions Wr{x,p) for some 
arbitrary chosen measurement numbers r. The evolution from a Poisson distribution in the 
number basis to the phase space distribution for the corresponding number state with n = 3 
is clear. As the phonon number converges, the phase uncertainty increases until we get full 
information about phonon number and phase uncertainty becomes complete. 



2. With mechanical damping 

Our simulations show a collapse to a mechanical phonon number state is possible even 
including loss from the second cavity. In an experiment the measurement induced collapse 
will be competing with thermal fluctuations and dissipation. The measurement induced 
collapse proceeds at a rate determined by the count rate which is bounded by k. We can 
estimate the rate of change of the phonon number variance due to the damping of the 
mechanics if we describe the mechanical dissipation using the usual weak damping master 
equation of quantum optics, 

^ = £p + 7^(7V + l)V[b]p^ + ^mNV[b'']pm (47) 

where the first term is the opto-mechanical part and 7^ is the mechanical damping rate and 
N is the mean thermal phonon excitation number at the mechanical resonance frequency. 
This equation indicates that phonons enter the mechanical resonator at a rate determined by 
7m iV and decay at a rate determined by ■jmiN + 1). The overall effect is a rate of increase 
for the number variance at a rate proportional to jmN. Thus the measurement induced 
state reduction will dominate the increase due to mechanical dissipation ii k » . 
We now consider the competition between thermalisation and the reduction in the number 
fluctuations due to successive measurements. 
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FIG. 9. Phonon number distribution function histograms after successive measurements for iV = 4, 
g = 1, K2 = 1, Ki = 0.2, 7 = 0.9 and different values of mechanical damping rate: (a) = 0, (b) 
7m = 0.00001, (c) 7m = 0.0001, (d) 7m = 0.001. In (b) we have the optimum value for mechanical 
damping rate for this scheme, so that we conditionally drive the mechanical resonator to a number 
state. We see the broadening effect of the mechanical damping rate as we increase 7m • 



If the system is started in the thermal equilibrium state that follows from Eq.(47), the 
initial state of the mechanical degree of freedom is a thermal state with the number distri- 
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bution 



1 °° N 

n=0 



The rate at with measurements can proceed is hmited by the time taken for each detection, so 
we will consider the case in which the single photon source is loaded immediately after each 
detection. We now include the effect of mechanical damping in the conditional Schroodinger 
for no detections up to time U, 

^ = _,(/^p _ pK^) + ^^(AT + l)V[b]p + ^mNV[b% (49) 



where K is given by the Eq.(23) and the damping superoperators are defined by Eq.d9 



Figure (|9j) shows the phonon number distributions after each readout for seventy mea- 
surements. When the mechanical damping rate is zero (figure 9(a)), the number distribution 
evolves from a thermal distribution into a n = 2 peak. If we increase the damping rate grad- 
ually, for K,2 ^ 10^^ we observe that the number distribution starts to broaden. It is clear 
from this discussion that the rate at which phonons enter and exit the mechanical resonator 
is the key limiting factor in reaching a Fock state. Thus the coupling of the mechanical 
resonator to its environment will need to be carefully engineered to make as small 

as possible, perhaps using an external phononic bandgap shield |10], or carefully engineered 
supports|32]- In the case of [7] a continuous-flow helium cryostat provides pre-cooling to 
achieve a bath occupancy of the a 3.68 GHz mechanical mode of iV < 100. 



III. CONCLUSION 



We have considered a quantum opto-mechanical system based on the coupling between 
a dual rail optical qubit code, formed of two optical cavity modes, and a single bosonic 
matter degree of freedom in the strong coupling limit for which the single photon opto- 
mechanical coupling rate g is of the order of the decay rate of the cavity k and shown 
how single photon detection can conditionally drive the mechanical resonator to a phonon 
number state provided the mechanical thermalisation rate is small enough. Despite 

the fact that in this realisation the interaction time between the mechanics and dual rail 
optical qubit is a random variable, we can found regimes in which the mechanical phonon 
number distributions become sharply peaked at a particular value of n. Successive photon 
counts, even though randomly distributed, provide a record of the interaction time in each 
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measurement which is sufficient to gain information about the state of the excitation number 
of the mechanics in a sequence of detection events. 

The achievement of the strong couphng regime will be a challenge. Several other exper- 
imental groups are working towards achieving strong coupling at the single photon level, 
[33H39] . The required values for the mechanical damping rate are achievable, for example, 
Chan et al. [10], have 7„ = 7.5kHz for k/2'k = 214MHz, g = l.lMHz and w„ = 5.1GHz. 
In the dimensionless units of this work these are, 7m//t = 3.5 x 10^^ for g/ n = 5.1 x 10~^ 
and Um/K = 30. The required value for the thermalisation rate 'jmN will be a challenge, et 
al. jlO] using a cryostat at 0.65K have 7miV/27r = 2MHz which with a optical decay rate of 
K = lOMHz corresponds to a dimensionless thermalisation rate of 0.3. We expect continued 
advances in the design and fabrication of electromechanical and optomechanical systems 
will enable sufficient isolation from the thermal environment to reach a regime where the 
measurement induced collapse of the phonon distribution can beat the broadening due to 
thermal effects. 
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